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Abstract 

The inclusive Ught hadronic decays of heavy quarkonia are studied within the framework of 
NRQCD at the leading order in v and up to the order of . With one-loop QCD corrections, the 
infrared divergences and Coulomb singularities in the decay amplitudes are proved to be absorbed 
by the renormalization of the matrix elements of corresponding NRQCD operators, and the infrared 
finite short-distance coefficients are obtained through the matching calculations. By taking the 
factorization scale to be 2mQ, the light hadronic decay widths are estimated to be about 274, 4.7, 
and 8.8 KeV for the r/c2, %2) and 77^2 respectively. Based on the above estimates, and using the 
El transition width and dipion transition width for the r]c2 estimated elsewhere, we get the total 
width of r]c2 to be about 660-810 KeV, and the branching ratio of the El transition r]c2 ^ he to 
be about (44-54)%, which will be useful in searching for this missing charmonium state through, 
e.g., the process ryc2 ^ ^ he followed by he ^rje- 
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I. INTRODUCTION 



The studies of production and decay mechanisms for heavy quarkonia provide important 
information on both perturbative and nonperturbative QCD. Based on the nonrelativistic 
(NR) nature of heavy quarkonium systems, an effective field theory, the nonrelativistic 
QCD (NRQCD) factorization formulism was proposed by Bodwin, Braaten and Lepage in 
1990s[l]. Within this framework, the inclusive decay and production of heavy quarkonium 
can be factorized into two parts, the short distance coefficients and the long distance matrix 
elements. Differing from the color-singlet model (CSM) j^, in the NRQCD factorization 
formalism, the heavy quark and antiquark pair annihilated or produced at short distances 
can be in both the color-singlet and the color-octet states with the same or different angular 
momentum quantum numbers l|, and the latter is known as the color-octet mechanism 
(COM). This mechanism has been used to remove the infrared divergences in inclusive P- 
wave charmonium production \^ and decay , 0] to give the infrared safe and model 

independent predictions. 

Recently, the inclusive light hadronic decays of ^Dj charmonium states were also studied 



within the framework of NRQCD factorization up to order a^|7|, l8|. The infrared divergence 

qI is removed by absorbing it into the matrix elements of the 



found in the CSM calculation 
color-octet ^Pj operators. Furthermore, the new contributions at order from the color- 
octet ^Pj and ^Si matrix elements enhance the decay widths of ^Dj states, and the numerical 
results are larger than those estimated in the CSM by several times in magnitude Q]. One 
can expect that a similar case will emerge in the inclusive light hadronic decay of 
charmonium, namely, the ric2 state. The difference between the ?7c2 and ^Dj states is that 
there are no infrared divergences in the inclusive decay width of ric2 in the CSM up to order 

[l^, and the numerical result is about 110 KeV 11 1. However, the infrared divergence 
will emerge again in the decay width of 7]c2 in the CSM at order a^, which needs to be 
removed by invoking the color-octet mechanism, i.e. by absorbing it into the corresponding 
color-octet matrix elements. 

On the other hand, the estimation of the inclusive light hadronic decay width of t]c2 is 
also important phenomenologically for probing this missing charmonium state. Quark model 



predicts its mass within the range 3.80-3.84 GeV 



12 



131 ]. which lies between the DD and 



the D*D thresholds. However, its odd parity {J^'" = 2 ^) forbids the decay to DD. As 
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a result, it should be a narrow state, and its main decay modes are the electric as well as 
hadronic transitions to lower-lying charmonium states and the inclusive light hadronic decay. 
Therefore, the study for the inclusive light hadronic decay of 77^2 in NRQCD factorization 
will provide important information on searching for this state in high-energy pp collision [ij] , 
in B decays 15|], in higher charmonium transitions, and in the low-energy pp reaction in 



PANDA at FAIR 



16| and in e+e" process in BESIII at BEPC [17 1. 



In this paper, we study the one-loop QCD corrections to light hadronic decay of 
within the framework of NRQCD factorization. The paper is organized as follows: after an 
introduction of the NRQCD factorization formalism in Sec. II, we calculate the decay widths 
up to 0{a^) in perturbative QCD in Sec. Ill, where both the real and virtual corrections 
are considered. Then perturbative NRQCD is applied to obtain the imaginary parts of the 
forward scattering amplitudes in Sec. IV. Combined with the QCD results, the infrared 
divergences are either canceled or absorbed into the long distance NRQCD matrix elements, 
and the finite short distance coefficients are obtained. Together with the long distance 
matrix elements estimated by solving the operator evolution equations, the decay width is 
determined. The numerical results and phenomenological discussions are given in Sec. V. In 
the last, we will give a brief summary for our results in Sec. VI. 



II. GENERAL FORMULAS 



There are four important scales in the heavy quarkonium system: the heavy quark mass 
rriQ, the typical momentum of the heavy quark or the inverse of the size of the bound state 
rnqv the binding energy rnqv"^ and the QCD scale Aqcd, while the dynamical property 
of the bound state is mainly determined by the latter three scales. Thus one can choose a 
cutoff fi\ with condition mq > /iA ^ rnqv {rngv'^, Aqcd) to integrate out the hard scale 
mg. Expanding the nonlocal effective action in power of v and writing the result in the two- 
component Pauli spinor space, one then can get the effective Lagrangian for NRQCD [l|: 

^NRQCD = flight + ^heavy + (l) 

Here, v denotes the relative velocity of the heavy quark pair in the meson frame. The average value of 
is about 0.3 for charmonium and about 0.1 for bottomonium ^|]. 
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where the Lagrangian Chght describes gluon and hght quarks. At leading order in v, the 
heavy quark and antiquark are described by Cheavy- 

Cheavy = i^\lDt + + x\lDt " ^)X, (2) 

where denotes the Pauh spinor field that annihilates a heavy quark, x denotes the Pauli 
spinor field that creates a heavy antiquark, and Dt and D are the time and space components 
of the gauge- covariant derivative D^, respectively. The relativistic corrections to Cheavy are 
included in the term 5C The most important correction terms for heavy quarkonium energy 
splitting are the bilinear ones: 



C2 



+ — V[V^n^D X - ^E X iD) ■ cr-^ + x"^(iD x - ^E x ?D) ■ ax] 

+ T^i^HaB ■ - xKgB ■ (3) 

where E' = G^' and B' = \e^^^&^ are the electric and magnetic components of the gluon 
field-strength tensor G^^ , respectively. 

In the Lagrangian Cnrqcd in dl]), there are still three low-energy sales: the soft scale 
mg V, the ultrasoft scale rriQ v"^ and the QCD scale Aqcd- The existence of multi-scales makes 
the power counting rules of NRQCD (the velocity scaling rules can not be homogeneous 
generally. More seriously, if one wants to do the NRQCD loop calculations in dimensional 
regularization scheme with Cnrqcd defined in ([1]), one will find that the hard scale can not 



decouple from the loop integrals and the power counting rules are violated^ inevitably 18|. 
These problems can be solved simultaneously by the method of regions [19,], which will be 
explained and applied in our calculations in Sec. IV. 

To reproduce the annihilation contribution to a low-energy QQ — > QQ scattering am- 
plitude in NRQCD, local four-fermion operators in 6C are needed, which have the general 
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form 



5Ca 

— fermion 

— rrin 



where C„ denotes regularized local four-fermion operators, such as ip^xx^i': and dn is the 
naive scaling dimension of the operator. The dependence on cutoff /^a of the operator On 



is canceled by that of scaleless coefficient /n(/iA), which can be computed by matching the 
full QCD onto the NRQCD as perturbation series in a^. 

In NR theory, the width of heavy quarkonium H is —2 times the imaginary part of the 
energy of the state, thus one has [l| 

T{H ^ LH) = 2lu,{H\5U.f,rrmon\H) = T^Ml^ (^H\On{fiA)\H) , (5) 

where LH represents all possible light hadronic final states, and the operator On here and 
afterward only denotes the one relevant to the strong annihilation of QQ. The NR normal- 
ization has been applied for the state \H) in ([5]). 

In order to calculate the coefficients of four-fermion operators in ([5]), the equivalence of 
full QCD and NRQCD at long distance is exploited. Since in construction, the coefficient 
fn is of short- distance nature and is independent on the long distance asymptotic state, one 
can get it by replacing the state \H) by the on-shell heavy quark pair state \QQ) with small 
relative momentum and matching the forward scattering amplitude of QQ QQ in full 
QCD onto that of NRQCD perturbatively. The matching condition is written as [l| 

A{QQ^QQ) =E^%^<^^i^«('"A)igg) . (6) 

pert QCD ^ pert NRQCD 

n V 

Since we only need the imaginary parts of the coefficients, optical theorem can be used to 
simplify the matching calculations. 

The physical state can be expanded in powers of v in the Fock space: 

\'D,) = 0{1)\QQCdP)) + Oiv)\QQCPl% + Oiv')\QQCsi;''^)) + 0{v'), (7) 

where the superindices [1] and [8] denote the color-singlet and color-octet, respectively. The 
contributions from the P-wave and S-wave Fock states to the annihilation rate of are at 
the same order of f ^ as that from the D-wave state, because their relevant operators scale 
and f ~^ relative to Oi(^D2)i as can be seen later. Other Fock states contribute at higher 
order of v"^. Therefore the light hadronic decay width of at leading order in v"^ can be 
described in NRQCD factorization framework as follows: 

r (1^2 ^ LH) = 

2 Im/(^Z.fV^^-'^^(^f-)'^^-^2Im/Cpf')^^^-'^-(^f^)'^^-) + 
rriQ ttIq 



where the four-fermion operators are 2l| : 

1 



OsCPi) 

01CD2) 



(9) 



where D = D - D and S'^ = {-^)^{D'D^ - lD^6'^). Since D^/ml scales as v'^, it can 
be ensured that the four terms in ([H]) are at the same order of v. 

The coefficients in can be obtained by applying the matching conditions (jS]) to ap- 
propriate QQ configurations. To subtract the full QCD amplitude of QQ state of particular 
angular momentum, the covariant projection method is adopted. In practice, the optical 
theorem can relate the imaginary part of the QCD amplitude ^ in ([6]) to the parton level 



2l| 



decay width [20|, 

T{QQ[n] ^ LFs) = ^{QQ[n]\0[n]\QQ[n])^^<'^Yl j l-^(WN - LFs)\^<l^, (10) 

where LFs denote the gluons or light quarks and [n] denotes the configuration of the QQ. 
The state IQQI'^]) has been normalized relativistically as one composite state with mass 
M = 2Eq, except that in the matrix element in ffTOl) . where the state is normalized non- 
relativistically to match the results in perturbative NRQCD conveniently. The super-index 
LO of the matrix element means that it is evaluated at tree level, and we always use the 
abbreviation {0[n])Lo to represent it in our calculations. Moreover, the summation/average 
of the color and polarization for the final/initial state has been implied by the symbol 

For spin-singlet states with L = 0, L = 1 and L = 2, the amplitudes defined in f fTOj) 
are given by |20| 



[1,8] 
So 



M{{QQ)[ 
MiiQQt^^ 




LFs) = a/— Tr[C[^'^]n°A^^ 



A^((gQ)[" 



D2 



LFs) 



iJD] 

2 




2 



M dqadqp 



(11) 



where denotes the parton-level amplitude amputated of the heavy quark spinors, 

and eir^ and e[^' are the polarization tensors for L = P,D states respectively. The factor 



2 _ V2M 



— . . comes from the normahzation of the composite state For 

color singlet and octet states, the color projectors are C'"*^! = and C^^l = y/2(Ta)ij 
respectively |20i]. The covariant spin-singlet projector 11° in flTT]) is defined by 

li' = Y,u{s)vm\,s;\Am- (12) 

ss 

The explicit form of 11° in D dimensions will be discussed in the latter subsection. 
The sums over polarization tensors for elT^ and e|^^ in D dimensions are: 



^S^^l^g* = 2^^'^a'^l3f5' + ^af5'^a'l3) — jj _ ^ ^a/S^a'fi'- (l^b) 

Here Haa' is defined as 

n, PaPa' /t a\ 

aa' = -Qaa' + ' (l^) 

where P is the total momentum of QQ., and = = 4£'q. 

Needless to say, the final result should be independent on the normalization convention 
of the QQ state. If one wants to apply NR normalization thoroughly in the calculations, one 



needs to eliminate the factors 1/(2M) in ( fTOl) and a/2/M in ( fTTl) . and then to replace the 
covariant spinors in (|T2l) with the NR ones with the normalization condition: u'^u = v'^v = 1. 



A. Discussions on 7^ scheme and projection operator 



We will do our calculations in dimensional regularization scheme both for QCD and 
NRQCD. Since we are only dealing with the spin-singlet Fock states, there will be the 
problem of definition of 7^ in D dimensions. In our calculation, the 't Hooft-Veltman (HV) 



scheme 



20 



22| is introduced: 



{7^7n = 0, /i = 0,l,2,3 
[75,^^] = 0, /x = 4,---,D-l. 



And the 7^ matrix can be represented as 23| : 
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4! 



e'^'^'"^7M7.7p7<.- 



(15) 



(16) 



The calculation involving 7^ is carried out in D dimensions, where e'^"''"' and 7^ are all 



defined in D dimensions. Other prescriptions may be found in literatures [241. |25|. I26l. 127 1. 

In four dimensions, the covariant spin-singlet projector 11° defined in (1121) can be given 
by (see, e.g., |28j) 

1 



(17) 



2x/2(Eq + mg) ' 2 

where q is half of the relative momentum of the heavy quark pair. The form in (I17p can not 
keep C parity conservation in D dimensions because [f + M)7^ can not keep an invariant 
form under charge conjugation transformation in D 7^ 4 dimensions in the HV scheme, 
which can be easily seen by applying (fT5l) . This problem can be solved by replacing it by 
the following two operators. For spin singlet states the spin projectors of incoming heavy 
quark pairs at any order in f ^ are given by 

1 ,t . Mfi + M)-i'> + -i\-f + M)] 



from 



2V2{Eq + mg) ' 2 
29|] and 

1 



2M 



^ - ttiq) (18) 



2v^(Eq + mg) ' 2 



+ M)^\-f + M) f 
2M2 ^ 2 



^q) (19) 



from 



30l |. The above two projection operators both give correct results and keep C parity 



conservation. 



III. FULL QCD CALCULATION 

In this section, we calculate the imaginary part of QQ forward scattering amplitude. 



or equivalent 
FeynArts 



y, the parton- level decay width F defined in (llOj) . In the calculation, we use 
31| to generate the Feynman diagrams and amplitudes and FeynCalc [32] for 
the tensor reduction. We regularize the ultraviolet (UV) and infrared(IR) divergence in 
dimensional regularization scheme and extend the covariant projection method into D = 
4 — 2e dimensions as has been mentioned. 

The leading order subprocesses in are the annihilations of QQ[n] into two gluons, 
where n can be any configurations of the Fock states listed in ([7]). The Feynman diagrams 
at LO of as are shown in Fig. 1. And the results in D dimensions are 
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rsTRRRRnr fcl 




FIG. 1: Feynman diagrams for ^Lj'^' ^ gg 



Born ' 



3 — 2e 



4e2 - I6e + 15 



where Bp = ^jj^ = ^ and $(2) is the two-body phase space in D dimensions: 
^ r/^o~o\ (-^Y- The first three resuhs in fl20l) are consistent with those in Ref. 

Stt i (2 — 2e) \ ttIq ' 



20 



2l|. 



At the Born level, there are no IR divergences in the results since both the two gluons 
should be hard in the rest frame of QQ. 



A. Real Corrections 



The real corrections to Born level subprocesses include the decays into ggg and qqg final 
states. The corresponding Feynman diagrams are shown in Fig. 2 and Fig. 3. For simplicity, 
unphysical polarization summation is used for final state gluons, so diagrams with ghosts in 
the final states must be included in calculation when three gluon vertex appears, in order 
to cancel the non-physical degrees of freedom to keep the full results gauge invariant. 



1- iQQ)„ 



[1,8] 



999 



Our results of S-wave configurations agree with those in 20|, |2l|] and are listed below: 



ggg) 



7T 
TT 



T ^orn{' S'i^ ^ gg)aM''){- + ^ + 



1 11 181 



23 



Born 



{'st^^gg)UM-'){- 



6e 18 24 
7 104 



-TT 



3e ^ 9 



— TT 



(21) 
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FIG. 2: Feynman diagrams for ^Lj'^^ ^ (7(jr(7 




FIG. 3: Feynman diagrams for Lj' qqg 



where /^(M^) = (i^)T(l + e). The D-dimension P and D -wave resuhs are: 



M2 

ri p[8] 



TV e"^ 27e lo2 



Both soft and colhnear IR divergences are there in the results in fl2T|) and fl22|) . and the 
square pole comes from the overlap of the soft and the colhnear regions. 



2. (QQ)i,[i,8i qqg 

Another subprocess of light hadronic decay is to qqg final states, and only two graphs 
make contribution to this subprocess (shown in Fig. 3). 
We get the following results: 



r 


(^4^1 - 


qqg) 


= NfTBorn 


. ^0 




-Tf{ 


2 


-) 
9 ^' 




r 


(^4^1 - 


ma) 


= A^/FBorn 


'1 q[8] 


7T K 


-Tf{ 


2 


16 
9 




r( 




ma) 


= NfTBorn 


'lpj8] _ 


71 K 


-Tf( 


2 
3e 


16 
9^' 




r( 


^D? - 


ma) 


= NfT Born 




aa) 




. 2 
^ 3e 


-) 

9 ^' 


(23) 
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where Nf is the number of hght flavor quarks. Nj = 3 and 4 for charmonium and bottomo- 

i, = r(l + e)r(l-e) c^l + e'^ 



nium respectively. Tp = ^ K = r(l + e)r(l — e) 1 + and the S-wave results agree 



with 



20 



2l|. 



There are only single poles of e in the results in ( l23i) and they can be identified as collinear 
ones. The absence of the soft IR divergence can be seen from the diagrams in Fig. 3. When 
the momentum of the real gluon goes to zero, it will decouple from the quark line as an 



eikonal factor 



2l|, then the results will be zero since QQ in spin-singlet can not couple to 



one virtual gluon. 

As will be seen later, the collinear divergences and partial soft IR ones in ( 12T1) . ( l22i) and 
(1231) are canceled by the virtual corrections to the Born level decay width. The remaining 
soft IR divergences are those in r(^P|^^ ggg) and T{^D2^ 999) from the first dia- 
gram in Fig. 2, which will be absorbed in the renormalization of the operators Oi^s^Sq) 
and Osi^Pi) in perturbative NRQCD. These are just the general results of the so-called 
topological factorization discussed in 



B. Virtual Corrections 



There are 23 virtual correction diagrams, including counter-term diagrams, divided into 
9 groups. Representative Feynman diagrams of each class are shown in Fig. 4. And the 
others can be found through reversing the arrows on the quark lines or exchanging the 
final state gluons. UV divergences are removed by renormalization. The definitions of the 
renormalization constant of QCD gauge coupling constant Qs = \/47raJ, heavy quark mass 
rriQ, heavy quark field ipQ, light quark field ipq and gluon field are: 

^° = Zg9s, rriQ = ZraQ-inq, ^'q = v^^^Q, = \^q^q, ^1 = \^Z^Af,, (24) 

where the superscript labels bare quantities, and Zi = 1 + SZi. The renormalized constant 
Zg is defined by minimal-subtraction (MS) scheme, and the others by the on-mass-shell 
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al 





a-i 




TnnrffTrtns^ 



a7 



()■ • 


Tnrcnrircnr 




(]' 






0, 





rtrsiyinnnr- 




a9 



FIG. 4: One-loop Feynman diagrams for {QQ)i^[i.8] — > gg 



(OS) scheme, similar to that in 



. Then the results are: 



5ZOS 



47r euv 
4:71 euv 

An euv 



- + 61n(2)+4), 
e 

h 

e 



sz. 



OS 
rriQ 



bZf' 



47r euv 

-6of^/.(M^)(— 
47r euv 



21n(2) + -) 



6 



TV; 
3 • 



(25) 



where bo 

We calculate diagrams one by one and summarize the results in the following form: 



r(i^[i,8] 



r(i^[i,8] 



n ^ — ' 



(26) 



where the results of T>k are listed in Table I-IV. We add the counter-term diagrams with 
the corresponding self-energy and vertex diagrams to show the explicit cancelation of the 
UV divergence. There are still IR divergences left, which will be canceled by those in the 
real corrections as we have mentioned. There are also the well-known Coulomb singularities, 
which have been regularized by the relative velocity v, in Table I-IV. These singularities can 
be absorbed by the corresponding matrix element through the matching condition ([5]). 
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TABLE I: Virtual corrections to (QQ)iam 99- 



Diag. 




ai+az 


/'I 1 1 1 it^ o^ 


ao+a4+ay 


2^ + ^( 2Gf &0+ 2 ) +''oln4^ + Cf( 8In2 4+^) 2 ( 4+ 


a5 


^^A(-^-i-2 + 21n2 + |7r2) 


a6 


^CA(^ + i + 2-41n2- 


a7 


<^^(i^ + I-2 + 21n2) 


a8 





TABLE IL Virtual corrections to {QQ)iqIs] — > 99- 


DiMg. 




al+a2 


CF(| + l + 61n2) 


a3+a4+a9 


-i^ + i(-2Cf-6o + ^)+6oln4^ + C^(-81n2-4 + ^)-%^(-4 + fJ) 


a5 


iCA(-^-i-2 + 21n2 + §7r2) 


mG 





a7 


(Cf -iCA)(g + i-2 + 21n2) 


a8 






C. Summary of the QCD results 

Combining the real and virtual correction results together and translating the parton- 
level decay width back to the imaginary part of the forward scattering amplitude, we get 



13 



TABLE III: Virtual corrections to iQQ)ipis] gg. 



Diag. 




ai+az 


Gf(- — o + iUmzj 




2^ £ i^Oi? + c»o + 2 j + "0 ill 4^ + 2 ioinzj ^ { b + 3 + zinij 


a5 


3Ca(-^- 7-5 + §7r2 + 51n2) 


a6 





a7 


{Cf-\Ca){^ + \-2 + 2\^2) 


a8 


Ca{-\ + '^) 



TABLE IV: Virtual corrections to {QQ)iy.[i] gg- 



Diag. 


V, 


al+a2 


Cf{\ - 10 + 22 In 2) 


a3+a4+a9 


- + 60 + Ca) + 60 In ^ + C^(10 + |7r2 - 38 In 2) - £i(-f + ^tt^ + 21n2) 


a5 


CA(-^-7-8 + |7r2 + 121n2) 


a6 




a7 


CF(g + i-2 + 4In2) 


a8 






the full QCD results up to C(a^): 



(2ImA(00|'S™] ^ QOCsfl])) = {3^(1 + ^) + ^[4(477 - 16^,) 

+12(33 - 2N,) In - 937r^]}(0('4'l)>io, 



(27a) 



lc[8]l 



pert QCD 



hual nag. 



12^; 



+ 



Set. 



432m| 



[16(153 - ANf) 



+12(33 - 27V^) In - 1297r^]}(0(^5f ^))^o, 



(27b) 
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a^[2(3(-8iV; - 211n(2) - 229) + Udn^) - 3{6Nf - 61)ln(^)] 



108m^ 



(27c) 



a3[4(128iV/ + 1008 ln(2) - 19509) + 192(A^/ - 9) ln(/V) + 72637r2] 



1620m^ 



(27d) 



where the states |gg[?^]) in the Lh.s. and the r.h.s. should be understood to have been 
normahzed under the same condition. Moreover, the equahties in (l27j) are independent on 
the normahzation conventions. That is, the state |gg[ra]) can be normahzed either relativis- 
tically or non-relativistically, either as a composite state or as a discrete state. Therefore it 
is convenient to use Im^ to do the matching calculations, and we will use the abbreviation 
Im^(n) to represent the amplitudes in fl27|) . 

The remaining infrared divergences and Coulomb singularities in (1271) will be precisely 
repeated in the radiative corrections of the matrix elements in perturbative NRQCD in next 
section. The finite short distance coefficients will be obtained after matching calculations. 



IV. NRQCD RESULTS AND OPERATOR EVOLUTION EQUATIONS 

In this section, we calculate the NRQCD corrections to the four fermion operators in D 
dimensions. As we have mentioned, we will adopt the method of regions 19| to avoid the 
mismatch of the loop momenta in different regions. Furthermore, each loop integral in this 
method contributes only to a single power in v, thus one can do the power counting before 
the integral has been explicitly done. 

Since the ultrasoft and the QCD scale are comparable, mgf^ ~ Aqcd, for both charmo- 
nium and bottomonium, there are only two low-energy scales to be considered in NRQCD, 
which satisfy the inequality rngv ^ rnqv'^. Thus, the nontrivial contributions to the 
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{T,p) i k 3 

soft : Qs ■■ ' = rh, A : ^AMAVAM/ = ^ 



+if. 



iT,P) 

potential : Q,, : » 



4 „ • — -i 4 . _ 'i^t? 



ultrasoft : Au : 'WWW = j^^^. 



FIG. 5: NRQCD Feynman rules for heavy quark and gluon propagators in different regions 



-— p, 



(a) (6) 
FIG. 6: NRQCD Feynman rules for heavy quark and gluon vertices 

NRQCD loop integrals come only from the following three regions: 

soft : : fco ~ 1^1 ~ ^qv, \E's : 7" ~ |p| ~ rnqv, 

potential : : /cq ~ rnQv"^, \k\ ~ mqv, \E'p : T ~ rnqv'^, \p\ ~ rnqv, (28) 

ultrasoft : : /cq ~ ~ niQv'^, 

where k,y and are the momenta of gluon field and heavy quark field respectively and 
T = po — rriQ. The loop momenta running in these regions scale as those of the corresponding 
gluons in dimensional regularization scheme. One can check that the other regions, such as 
that with ko ~ rnqv and k ~ rngv'^, have no contributions to the NRQCD loop integrals 
in dimensional regularization scheme. In the view point of effective field theory, the five 
modes defined in ( 128|) can be all treated as the effective fields in NRQCD. Only after parting 
these low-energy modes sufficiently like what has been done in ( 128|1 . the homogeneous power 
counting rules can be gotten. 



In practice, we use the NRQCD Feynman rules 



34| derived in Coulomb gauge for the 



three regions (or the five low-energy modes) in our calculations. These rules are shown in 
Fig. 5 and 6, where = 5^^ — The Feynman rules for anti- heavy quark could be 

obtained by charge conjugation symmetry. 
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FIG. 7: NRQCD Feynman diagrams for LO Matrix Elements 



(T,p) (T',p') (T,pj {T\p') 




FIG. 8: NRQCD Feynman diagrams for NLO Matrix Elements 



The Coulomb singularities calculated in full QCD theory correspond to the potential 
region, while the soft divergences to the soft one. The LO Feynman diagrams for matrix 
elements are shown in Fig. 7. The onshell external quark lines lie in potential region. At 
NLO in as, only six classes of Feynman diagrams shown in Fig. 8 need to be calculated for 
our purpose 1,81]. The first four diagrams (a)-(d) have inner gluon lines connecting with one 
mcomiug quark liue and one^outgoiug quark hue. aud the soft regiou will grve the lowest 
order nontrivial result in v [8]. In the last two diagrams (e) and (f) the inner gluon line 
joints two incoming or outgoing quark lines, and only the potential region has non-vanishing 
real value |8]. The self-energy diagrams in the external legs are dropped in accordance with 
the on-shell renormalization scheme used in the full QCD calculation. 

We present here the detailed calculation of the NLO correction to the P-wave octet 
operator O(ipfl). The LO result (C(^Pf ^))Born is trivial. Using the Feynman rules for 
propagators of heavy quark and gluon in soft region and for the heavy quark gluon vertex 
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between potential and soft regions, the loop integral of diagram (a) is: 
^ igl f d'^k p-p'-(p-k)(p'-k)/k2 1 1 



(29) 



tUq J (27r)-^ kl — 'k^ + ie ko — ie ko — ie' 

After performing the contour- integration of A;o = |k| — ie, we get^ 

r _ g's f d^-'k p-p--(p-k)(p--k)/k^ 

which is both infrared and ultraviolet divergent. The integral in (150]) is scaleless, so vanishes 
in dimensional regularization. That is, the UV pole will be canceled by the IR one. But the 
result is nontrivial: 

' ^ -)P-P'. (31) 



Sttuiq euv e 

The integrals of (b)-(d) in Fig. 8 could be evaluated in the same way, and their results are: 

A-. = Tr^(— --)p-p. (32) 

Svrm^ tuv e 

Making use of the Feynman rules for heavy quarks and gluon in potential region, we 
obtain the loop integral of diagram (e): 

where T = Integrating fco? we have 

= ^^"^^ J (2^k^ k^ + 2p.k-.e - ^^'^ 



This integral could be performed directly by introducing v = and we get the Coulomb 
singularity: 



= ^(1 - -(- - In(^) - ie)). (35) 

The integral of diagram (f) gives the same Coulomb singularity but with opposite sign in 
the imaginary part. 



b 



Since the quark propagator poles should be taken into account in the potential region, one only needs to 



evaluate the contribution from the gluon pole in the soft region IS | . 
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The color structures of diagrams (a,c), (b,d) and (e,f) are obtained by color decomposition 
and are listed below respectively: 

1 1 _ 2 

rpb^j.aj.b ^ ^j.a ^ (^^^ _ }lCa)V2T^ ® V2T^. (36) 

Summing over each integral multiplied by the according color factor we get NRQCD 
matrix element of P-wave operator at NLO, which is UV divergent and needs to be renor- 
malized: 

(O0CPfl))^LO = {(l + ^(CF-|C^))V2T^®v^T^+ ^ ^ 

An (-1^)^'' _ (37) 

where we have used (C°(^Pf ^)) = ((!)°(^Pi)) ® and the superscript means the 

bare operator. Before doing the operator renormalization, we first re-express the bare result 
as 

m'P?^))r,Lo = (1 + ^{Cf - |C^))(0(^Pf' ))lo 

From the above equation we can see that the color-octet P-wave operator is mixed with the 
color-singlet D-wave operator at NLO in a^. We define the renormalized operator 0^{^pI^^) 



through 



33| 



{OrPl'^))NLO = {OXP?)hLO + '-^t^iT^ + 1^471 - 7e)(C^(0(^/^?1))lo 

+BF{OeDf^))Lo). 



Here the MS renormalization scheme is adopted. The matrix element of the renormalized 
operator is UV finite, but still has IR divergence term which will cancel the infrared divergent 
D-wave full QCD result. And it also has the Coulomb singularity, which is the same as that 
appearing in the full QCD virtual correction: 

{OXpI'^))nlo = (1 + ^{C, - ICa))WpI'^))lo 
+^^i-7 - ln47r + 7e)(CHO(1dW))z.o + B,{oM^)),o). 

Here, the matrix element of (9(^1^2^') is at higher order in v in our case and therefore can 
be eliminated. The matrix elements of S-wave singlet and octet operators and that of the 
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D-wave singlet operator could be computed in the same way: 



_i_ln4vr + 7E)WPfl))LO, 



1C'[1]^ 



I 1 f 1 

~'~27Vc 37rm2 



-l-ln47r + 7i.)(0(^Pf^)LO + ..., 



37rm^ 



2t; 



(41) 

(42) 
(43) 



where denotes terms at higher order in v. 

Finally combining the matrix elements given above with the short distance coefficients 
accordingly, we get the forward scattering amplitudes for ^Lj states computed by NRQCD 
effective theory, which are summarized below: 



(2W(i5W)) 



21111/(15™: 



pert NRQCD 



2v 



LOi 



(44a) 



{2\mA{'sf)) 



2Im/(i5, 



■[l + {C,-^CA)^]{0{'sf))Lo 



pert NRQCD 



2v 



(44b) 



(2W(ipfl )) 

pert NRQCD 

1 4a, 2Im/(i5'f'; 
'2Nc37rm^Q e 



AasBF2lmfi^Sj^^] 



(44c) 



{2lmACD?)) 



_ 2Im/(iZ}M; 

pert NRQCD mg 

lp[8]N 



2t; 



AasCp2lmf{'Pf 



(44d) 



37rm^ e 

Setting expressions in equal to those in fl27p respectively, and expanding Im/„ in 
power of as, we obtain the IR finite short distance coefficients up to 0{al): 

,2 



2lmfCS, 



/ 



+ -^(4(477 - 16Nf) + 12(33 - 2Nf) ln( 



Ami 



937r' 



(45a) 



2Im/(i5, 



lQ[8h 
/ 



^ + ^(16(153 -4iV,) + 12(33 -2iV,)l„(£^ 



1297r^ 



(45b) 
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2Im/(ipfl) = ^ + ^{-3(6iV^ - 61) ln(£^) + 2[-{2ANj + 631n(2) + 725) 

+ 1197r2 + 1141n(— )]}, (45c) 

1 R 2 3 

2Im/(^DW) = + [78720 - 512iV/ - 72637r2 - 4032 ln(2) + 2880 ln(-^) 

-192(iV;-9)ln(£^)], (45d) 

where the short distance coefficients of P-wave and D-wave are /iA-dependent. Their /za- 
dependence will be canceled by that of the corresponding renormalized operators, which 
could be obtained by finding the derivative of both sides of ( I40|41|l42p of /xa. For Born 
quantities, — = 0. Then we obtain the renormalization group equations at leading 

order in v and a,: 



(iln/iA S^'^^ 



d In /iA 2Nc 3vrm^ 

^ (46) 

The solutions of the matrix elements {1^D2\0{^a)\1^D2) in heavy quarkonium D-wave state 
are: 

(l^D^IO^l^Pf )(/iA)|l^/^2) = -^ln^i^(liD,|0(iDW)|liD,), 
(l^A|O^C5W)(M|liZ^,) = ^(:^ln^4^)^(l^D,|OrDW)|l^D,), 
(l^D,|O^C5r)(;.A)|l^/^2) = ^(;r4l-ln^4^)^(l^/^2|OC/^™)|l^/^2), 

(47) 

where the initial matrix elements like (1^Z^2|C^'^(^-Pi^')(/^Ao)I-'^^-^2) at /xaq = ttT'QV are elimi- 
nated at LO in V ^. 



V. NUMERICAL RESULTS AND PHENOMENOLOGICAL DISCUSSIONS 



The long distance matrix element of D-wave four-fermion color-singlet operator in the 
hadron state is related with the second derivative of radial wave function at the origin 
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FIG. 9: Renormalization scale dependence of the decay width of charmonium state 1^D2 to LH 



through the following relation: 



{n'D,\0{n'D,)\n'D,) = '^'"^g^^^^' = ml^Hnn ■ (48) 

ln[ll^ 



Combining the leading order coefficient 2Im/(^D2 )lo = 167ra^/45 given in ( 145d[) and the 



color-singlet matrix element given in fHHj) . one can reproduce the decay width of state 
in the CSM at leading order in |lo| |: 

row("'C.-M) = ^i^H5)^. (49) 

However, there are contributions from the color-octet Fock states in (JTj) in NRQCD even at 
the order of af. The matrix elements of the P-wave octet operator and S-wave singlet as 
well as octet operators in the bound state could be estimated through the solutions of 
the operator evaluation equations in (H7I) . 

The region of validity of the evolution equation is chosen as follows: the lower limit 
A^Ao = ^QtJ and the upper limit fiA of order itlq. For convenience, we take the factorization 
scale to be the same as the renormalization scale /i of order rriQ. We choose the pole mass 
TTir = l.SGeV, = 0.3, fiAo = rricV, ua = 2mc, as{2mc) = 0.249, Nf = 3, Aqcd = 390MeV, 



|2 



Hni = = 0.786 X lO-^GeV |35| for charmonium, and nib = 4.6GeV, v"^ = 0.1, /xao 

mhV,fXA = 2m,,asi2m,) = 0.180, Nj = A, Aqcd = 340MeV, Hoi = '"^^f ' = ^f^^ ^ 
10-^GeV for ID states and Hd2 = ^^'f^'g^f = 0.750 x lO'^GeV for 2D states [331 for 
bottomonium. The /x dependence curves of the decay widths are shown in Fig. 9 and Fig. 
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FIG. 10: Renormalization scale dependence of the decay width of bottomonium state 1^D2 to LH 

TABLE V: Subprocess decay rates of 1^D2 charmonium, where nic = 1.5GeV, v"^ = 0.3, fiAo = rricV, 
= 2mc and as(2mc) = 0.249. 



Subprocess 


LO(KeV) 


NLO(KeV) 


CD2)i ^ LH 


54.7 


75.1 


(iPi)8 ^ LH 


66.6 


132 


(i5o)8 ^ LH 


15.0 


31.3 


(iSo)i ^ Lii- 


19.2 


36.1 



10. When n — 2mc for cc systems and 2mf, for bb systems, we get the predictions at 0{al): 

rc(l^£>2 ^ LH) = 274KeV, 
rB(l^D2 ^ i-^) = 4.70KeV, 

rB(2iL>2 ^ Li/) = 8.78KeV. (50) 

The LO decay widths for charmonium and bottomonium 1^D2 states are 155 KeV and 3.22 
KeV. Therefore the NLO QCD corrections contribute enhancement of factor 1.8 and 1.5, 
respectively. 

For the 1^D2 charmonium state r)c2, the numerical values for all subprocesses are also 
hsted in Table V. One can see from this table that the contributions from the Fock states 
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other than |^-D2^') are dominant in the decay width, and the total result is about 1-3 times 
larger than that in CSM even at leading order in a^. 

For the phenomenological analysis of ric2, we vary the renormalization/factorization scale 
from 2mc to rric and get F (77^^ LH) = 274-392 KeV. The electric transition rate r{?7c2 — *■ 



7/ic) = 339-375 KeV [l3l] and the dipion transition rate T{'ric2 tIcTctt) ~ 45 KeV have 
been estimated elsewhere in the literature. 

As emphasized before, the r]c2 should be a narrow state, since its mass and quantum 
numbers forbid it to decay into charmed meson pairs DD and D*D. Therefore, the main 
decays modes of ?7c2 are expected to be the electric as well as hadronic transitions to lower- 
lying charmonium states and the inclusive light hadronic decay. With all these decay widths 
given above, we get the total width of ?7c2 to be about 660-810 KeV, and the branching ratio 
of the electric transition to be 

S(r7e2^7M = (44-54)%, (51) 

which provides important information on probing this missing state. In practice, one can 
search for ?7c2 through the cascade decay ?7c2 ih^ — > 7777^ — > j'jKKtt with branching 



ratios B{hc "yrjc) ~ 0.4 [ll|, |37j and B{ric KKir) ^ 7% [36]. Similar decay chains can 
also be used to search for the 77^2 • 

The production rates of ric2 are expected to be generally low in many processes, because 
the rates are suppressed by the small values of the second derivative of the wave function 
at the origin, and also by its spin-singlet nature, which forbids r]c2 to couple to a photon, 
or to be detected from the El transitions of higher spin-triplet charmonia. Nevertheless, 
efforts should be made to find this very unique missing charmonium state. Hopefully, the 
study for the inclusive light hadronic decay of ric2 in NRQCD will provide useful informa- 
tion on searching for this state in high-energy pp collision [l^, in B decays [l^, in higher 



charmonium transitions, in e~^e process in B 



low-energy pp reaction in PANDA at FAIR 16 1. 



vSIII at BEPC [l7|, and particularly in the 



VI. SUMMARY 

In this paper, we calculate the inclusive light hadronic decay width of the heavy 
quarkonium state up to order of and at the leading order in v within the framework of 
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NRQCD. We find that the inclusive decay widths into hght hadrons via gluons and hght 
quarks at order of in QCD suffer from both IR divergences and Coulomb singularities, but 
they can be absorbed into the renormalization of the matrix elements of the four-fermion op- 
erators in NRQCD precisely. Therefore, after matching the full QCD onto NRQCD, the IR 
divergent part is removed, and IR finite short-distant coefficients are obtained, and the de- 
pendence on the factorization scale of the coefficient is canceled by that of the corresponding 
matrix element with the renormalization group analysis. 

At leading order in a^, the result in the CSM can be reproduced but there are many 
other contributions, such as that from color-octet P-wave operators, which will enhance the 
width in CSM by several times in magnitude even at the leading order in a^. Furthermore, 
the NLO results give extra enhancement factors of 1.8 for ric2 and 1.5 for rih2 relative to the 
LO ones, respectively. By choosing the factorization scale as 2mQ, the light hadronic decay 
widths are found to be about 274, 4.7, and 8.8 KeV for the 77^2, Vb2, and 77^2 respectively. 
Based on these estimates, and using the El transition width and dipion transition width 
for the 7]c2 estimated elsewhere in the literature, we get the total width of r]c2 to be about 
660-810 KeV, and the branching ratio of the electric transition ric2 7 /ic to be about 
(44-54)%, which will be useful in searching for this missing charmonium state through, e.g., 
the process ?7c2 — > 7 /ic followed by he — ^ 777c- 
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